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Abstract 

O . In analogy to the Altarelli-Parisi equation for the quark and gluon helicity 

contributions to the nucleon spin, we derive an evolution equation for the 
quark and gluon orbital angular momenta. The solution of the combined 
^ ■ equations yields the asymptotic fractions of the nucleon spin carried by quarks 

and gluons: 3nj/(16 + 3nj) and 16/(16 + 3nj), respectively, where rij is 
the number of active quark flavors. These are identical to the well-known 
asymptotic partitions of the nucleon momentum between quark and gluon 
contributions. We show that the axial-anomaly contribution to the quark 
helicity is cancelled by a similar contribution to the quark orbital angular 
momentum, making the total quark contribution to the nucleon spin anomaly- 
free. 
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What are the composition of the nucleon spin in terms of its quark and gluon con- 
stituents? In the last few years, the answer has become the holy grail for the nuclear and 
particle spin community. The motivation is quite clear: given Quantum Chromodynamics 
(QCD) is difficult to solve, such information gives us valuable insights into the nonperturba- 
tive wavefunction of the nucleon. A satisfactory understanding of the spin structure would 
be to know, for instance, how much spin of the nucleon is carried, respectively, by the quark 
and gluon helicities and orbital angular momenta, in the same way as one now understands 
how the mass of the nucleon is partitioned among contributions from the quark and gluon 
kinetic energies, quark masses, and the trace anomaly (T|. 

The theoretical basis for separating the spin of the nucleon into different contributions 
begins with the QCD expression for the generators of Lorentz transformations f2|, 

jw = J d 3 xM^ u , (1) 

where M atJil/ is the angular momentum density, 

M afMU = i^ a {x^d u - x^^ifj + y>'f['f,'f]i> 

-F aa (x^d v - x v d^)A a - F afl A u + F au A» , (2) 

where color indices are implicit. According to definition, a nucleon moving in the z direction 
with momentum P^ helicity 1/2 satisfies, 

J 12 \P+) = \\P+) ■ (3) 
Thus one can write down a spin sum rule, 

I = (P+ |j 12 |p+)/(p+ |p+) 

= ^A£ + Ag + L q + L g , (4) 
where the matrix elements are defined as (7 s = 27°7 1 7 2 7 3 ), 

d 3 x V>7 3 7 5 V' \P+) , 



AS = 


(P + 


\$3q\ 


P+) = 


(P + 


Ag = 


(P + 




\P+) = 


(P + 


L q = 


(P + 


\^3q 


\P+) = 


(P + 


L 9 = 


(P + 


\ L 3g 


\P+} = 


(P + 



J d 3 x (E l A 2 - E 2 A 1 ) \P+) , 

J d 3 x i^^d 2 - x 2 d l )il) \P+) , 

J d 3 x E\x 2 d 1 - x^d^A 1 \P+) , (5) 



where for simplicity we have neglected the normalization of the state. It is clear from the 
above that AE and Ag are the quark and gluon helicity contributions to the nucleon spin, 
and L q and L g are the quark and gluon orbital angular momentum contributions. Apart 
from £>3 q , the other three operators £>3 g , L 3q and L 3g are not manifestly gauge invariant, and 
thus a decomposition of the nucleon spin is in general gauge-dependent. Furthermore, the 
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matrix elements depend on the choice of Lorentz frame. Only in light-front coordinates and 
light-front gauge |3) Ag is the gluon helicity measured in high-energy scattering processes. 
We henceforth work in this coordinates and gauge 0] (the index in Eq. (1) is now replaced 
by+). 

The individual operators in J 12 are not conserved charges, and hence their matrix el- 
ements are generally divergent. They can be renormalized in a scheme and the renormal- 
ization introduces a scale dependence. The scale dependence of AS and Ag obeys the 
well-known Altarelli-Parisi (AP) equation, which in the leading-log approximation is 0, 

(6) 

where t = lnQ 2 /AQ CD , Cp = 4/3 and @ = 11 — 2rif/3 with rif the number of quark flavors. 
The square matrix in the right-hand-side is called the splitting matrix. In this paper, we 
derive an equation for the leading-log evolution of the quark and gluon orbital angular 
momenta. The result is 







We call the first term in the right hand side the homogeneous term and the second the 
inhomogeneous term (with corresponding splitting matrices). The inhomogeneous term 
was first studied by Ratcliffe [0, from which he concluded that orbital angular momentum 
plays an essential role in the nucleon spin. At operator level, this means that the orbital 
angular momentum operators contain leading-twist contributions. However, our result for 
the inhomogeneous part of the L q — L g evolution disagrees with Rat cliff e's. Furthermore, 
the homogenous term we have is new. 

In the remainder of the paper, we sketch our derivation of Eq. (0) and find its general 
solution. At the end of the paper, we go beyond the leading-log approximation and discuss 
the "anomaly cancellation" in the angular momentum operator J 12 . To keep the physics 
clear, we will follow as closely as possible the original language of Altarelli and Parisi in light- 
front coordinates although the whole discussion can be made consistently in the language of 
operator mixing. When necessary, we supplement our discussion with the matrix elements 
of angular momentum operators in composite parton states. 

To begin, we review the standard derivation of the AP equation. Consider a parent quark 
with momentum^ = (p~ = 0,p + ,p± = 0) and helicity +1/2, splitting into a daughter gluon 
of momentum fc M = (k~ ,xp + ,k±) and a daughter quark with momentum (p — k)^. Only 
3-momentum is conserved during the splitting. The total probability for the splitting is, 
dx(l + (1 — x) 2 )/x. [A multiplicative factor a s hiQ 2 / /z 2 is implied when we talk about 
probability, /i 2 here is an infrared cut-off which defines the nature of the parent quark. It 
must be large enough so that perturbative QCD is valid. Q 2 is an ultra-violet transverse 
momentum cut-off which defines the scale of daughter partons.] Since the quark helicity 
is conserved at the leading-log, we therefore have the item in the upper-left corner of 
the AP splitting matrix. The helicity of the daughter gluon can either be +1 or —1. The 
probabilities for both cases are Jq dxl/x and Jq dx(l — x) 2 /x, respectively. The gluon helicity 
produced in the splitting is just Jq 1 dx(l — (1 — x) 2 )/x = (3/2)Cp, which is the element in 
the lower-left corner of the AP splitting matrix. 
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A further consideration of the above process leads to an inhomogeneous source for the 
orbital angular momenta. Since the total angular momentum in the z direction is conserved 
in the splitting, orbital angular momentum has to be produced to cancel the helicity of the 
daughter gluon. This means that total orbital angular momentum carried by the daughter 
quark and gluon is — (3/2)C F . QCD determines partition of this between the quark and 
gluon according to the matrix elements, 

(p+\L 3q \p+), (p+\L 3g \p+), (8) 

where \p+) is the quark-gluon state produced in the splitting. To calculate them, we start 
with less singular off-forward matrix elements and consider their forward limit B . We find, 



(p+\L 3q \p+) = 


-C F jV [ 
Jo 


dxx { 

\x 




-\c F N , 




(p+\L 3g \p+) = 


—C F J\f [ 
Jo 


dx(l — x) 




-\c F N , 







1 (1-x) 



X X 



(9) 

where IV = 2p + (27r) 3 <5 3 (0)(a s /27r) lnQ 2 //z 2 . The coefficients in front of IV give the first 
column of the inhomogeneous splitting matrix in Eq. (|7|) . Interestingly, the sharing of the 
orbital angular momentum is done according to the x and 1 — x moments of the polarized 
gluon density in the parent quark. 

A similar discussion leads to the second columns of the AP splitting matrix and the 
inhomogeneous splitting matrix in Eq. (|7]). Let us remark here the physical origin for the 
well-known result that the gluon helicity increases logarithmically in the asymptotic limit 
jTj. When a gluon splits, there are four possible final states: 1) a quark with helicity 1/2 
and an antiquark with helicity —1/2; 2) a quark with helicity —1/2 and an antiquark with 
helicity 1/2; 3) a gluon with helicity +1 and another with helicity —1; 4) two gluons with 
helicity +1. In the first two processes, there is a loss of the gluon helicity with probability 
(rif/2) Jq 1 dx(x 2 + (1 — x) 2 ). In the third process, there is also a loss of the gluon helicity, 
with probability Jq dx(x 3 /(l — x) + (1 — x) 3 /x). In the last process, there is a gain of the 
gluon helicity with probability Jq dxl/(x(l — x)). When summed, the gluon helicity has a 
net gain with probability 11/2 — rif/3 = (3q/2 in the splitting. Thus at increasingly smaller 
distance scales, gluons split consecutively and the helicity builds up logarithmically. 

Now we turn to the homogenous term in the orbital angular momentum evolution equa- 
tion. To calculate the splitting matrix, we consider the matrix elements of L 3q and L 3g in a 
parton state with non-vanishing transverse momentum. In a bare parton state, we have, 

(p' + |L 3 Jj9+) = 2p + (27r) 3 ^j9;^ + ^^5 3 (p / -p) . (10) 

The derivative on the 5-function means that when the distribution is convoluted with a test 
function, the derivative will be taken of the test function. When calculating the matrix 
element in the composite parton states from the parton splitting, we have 
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ip' + \L Sq \p+) = 2p + (2rt) 3 <j ) (p\p) (-ip'i-^r + <*V ~ p) • (11) 

When the distribution is convoluted with a test function, the derivative will be taken of 
<f)(p',p) and then of the test function. The first term represents generation of the orbital 
angular momentum from the parton helicity discussed previously. The second term has the 
same structure as the basic matrix element in Eq. ([H^) and represents the self-generation of 
orbital angular momentum in the splitting. 

According to the above recipe, we calculate the gluon orbital angular momentum gener- 
ated from the quark orbital angular momentum, 



Z7T jjL z JO X 



(I 2 ) 



— In — C F - 

2ir /X 2 3 

Interestingly, this is the same as the fraction of the quark momentum carried by the daughter 
gluon in the infinite momentum frame. The net loss of the quark orbital angular momentum 
in the splitting must be the negative of the above. These numbers form the first column of 
the homogeneous splitting matrix. A similar calculation yields the second column. 

The solution of the evolution equation can be obtained straightforwardly. First, let us 
write down the well-known result for the quark and gluon helicities [§,0) 

AS(t) = const , 

. , , 4AS t / 4AS\ , , 

The second equation exhibits the famous behavior of the gluon helicity: increasing like In Q 2 
as Q 2 — > 00. The coefficient of the term depends on the special combination of initial quark 
and gluon helicities, which is likely positive at low-momentum scales according to the recent 
experimental data on AS M and theoretical estimates for Ag The solution for the 

orbital momenta is, 



3n 



f 



16 + 3n 



LJt) = --AE + - 3n/ + (t/to) - Wn,)/^) ( L (0) + i AE _ i 
qK 1 2 216 + 3n/ w ' \ qy ' 2 2 

Given a composition of the nucleon spin at some initial scale Qq, the above equations yield 
the spin composition at any other perturbative scales in the leading-log approximation. 
From the expression for L g (t), it is clear that that the large gluon helicity at large Q 2 is 
cancelled by an equally large, but negative, gluon orbital angular momentum. 
Neglecting the sub-leading terms at large Q 2 , we get, 

1 1 3n f 

J = L a + -AS 



2 2 16 + 3n/ 

1 16 

J - L + Aq = . (15) 

9 9 y 2 16 + 371/ K J 
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Thus partition of the nucleoli spin between quarks and gluons follows the well-known par- 
tition of the nucleon momentum fl l|] ! Mathematically, one can understand this from the 
expression for the QCD angular momentum density = T m x^ — T^x a . When 

and T Q/3 are each separated into gluon and quark contributions, the anomalous dimensions 
of the corresponding terms are the same because they have the same short distance behavior. 

It is interesting to speculate phenomenological consequences of this asymptotic partition 
of the nucleon spin. Assuming, as found in the case of the momentum sum rule [12], that 



the evolution in Q 2 is very slow, then the above partition may still be roughly correct at low 
momentum scales, say, Q 2 ~ 3 GeV 2 . If this is the case, from the experimentally measured 
AS we get an estimate of the quark orbital contribution at these scales, 

L q = 0.05 ~ 0.15 . (16) 

To find a separation of the gluon contribution into spin and orbit parts, we need to know Ag, 
which shall be measurable in the future |H| . However, if Q 2 variation is rapid, the asymptotic 



result implies nothing about the low Q spin structure of the nucleon. Unfortunately, no 
one knows yet how to measure L q to determine the role of Q 2 variation. 

So far we have considered the leading-log results. It is well-known that at the next- 
to-leading-log level, the matrix element of the axial current in a gluon state contains an 



anomaly [0,|14[ . Because of this anomaly, the axial current acquires an anomalous dimension 
starting at the two- loop level ^ JTHjTS] . Also because of the anomaly, it was suggested that 
the quark helicity for light flavors be modified to subtract the anomaly contribution |T7||7|,p!4| . 
Since the axial current appears in the conserved charge J 12 , the anomaly contribution to 
AE must be cancelled by matrix elements of some other operators in J 12 . Below, we show 
such cancellation comes from the quark orbital angular momentum L 3q . 

To investigate the anomaly cancellation, we consider a gluon moving in the z direction 
with momentum p^ and helicity +1. The gluon can virtually split into a quark- ant iquark 
pair. For zero quark masses, the quark and antiquark have opposite helicities due to chirality 
conservation. Thus the helicity of the gluon is entirely transferred to the orbital angular 
momentum of the pair. However, in the limit of the quark transverse momentum going to 
infinity, the chiral symmetry is broken by the anomaly. As a consequence, the quark and 
antiquark can have same helicity. Indeed, the calculation of Carlitz, Collins and Mueller 
shows []14| , 



1 - a 
(P + \^1 + 7^\P+) = ~ x 2p + . (17) 

The sign indicates that the quark and antiquark prefer to have helicity —1/2. To preserve 
the total angular momentum, the orbital motion of the pair has the same probability to 
carry two units of angular momentum. Indeed, an explicit calculation yields, 

(p + \L 3q \p+) -(p+ \p+) = x 2p+(27r) 3 5 3 (0) . (18) 

According to the above, if one redefines the quark helicity by subtracting the anomalous 
gluon contribution |7,14], one must also redefine the quark orbital angular momentum by 



adding to it the same amount. Thus if a large Ag helps to restore the measured AS to the 
quark model result, the same effect will make L q in Eq. (16) negative. However, the fraction 
of the nucleon spin carried by quarks is invariant under such redefinition. 
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To summarize, we have derived an evolution equation for the quark and gluon orbital 
angular momenta in QCD. In the asymptotic limit, the solution of the equation indicates that 
the quark and gluon contributions to the nucleon spin are the same as their contributions 
to the nucleon momentum. Furthermore, both the gluon orbital and helicity contributions 
grow logarithmically at large Q 2 and with opposite sign, so their sum stays finite. We show 
that the anomaly in the quark helicity is cancelled by a contribution in the quark orbital 
angular momentum, so the total quark contribution to the nucleon spin is anomaly-free. We 
also discussed phenomenological implications of our results, in particular, about the size of 
the quark orbital angular momentum. 
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